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Receptivity of Supersonic Boundary Layers Due To Acoustic 
Disturbances Over Blunt Cones 
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NASA Langley Research Center, Hampton, VA 23681 

Receptivity and stability of supersonic boundary layers over a 5-degree straight 
cone with a blunt tip are numerically investigated at a free stream Mach number of 
3.5 and at a high Reynolds number of 10 6 /inch. Both the steady and unsteady 
solutions are obtained by solving the full Navier-Stokes equations using the 5 th -order 
accurate weighted essentially non-oscillatory (WENO) scheme for space 
discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta 
scheme for time integration. The linear stability results showed that bluntness has 
less stabilizing effects on the stability of boundary layers over cones than on flat 
plates and wedges. The unsteady simulations of the interaction of plane three- 
dimensional acoustic waves with the cone showed that the modulation of wavelength 
and the generation of instability waves first occurred near the leading edge in the 
plane where the constant acoustic phase lines are perpendicular to the cone axis. 
Further downstream, this instability region spreads in the azimuthal direction from 
this plane. 


Introduction 

This paper is concerned with the receptivity and the stability of supersonic boundary layers over blunt 
axi-symmetric cones at zero angles of attack. It is now understood that transition from a laminar to a 
turbulent state originates from the internalization of external disturbances into the boundary layer 
(receptivity) and the exponential growth (linear instability) of these internalized disturbances inside the 
boundary layer. The second process is governed by the linearized stability equations. The stability 
characteristics of compressible boundary layers have been thoroughly investigated 1 . The supersonic 
boundary layers are unstable to three-dimensional disturbances called the first mode. The phase velocity 
of the unstable first mode disturbances increases from LL - a^ near the lower neutral point to C s near the 
upper neutral point, where LL is the free stream velocity, a. x is the free stream acoustic velocity and C s is 
the mean velocity at which the generalized inflection point occurs. The unstable region and the 
amplification rate of the disturbances depend on the flow parameters, Mach number and unit Reynolds 
number, on the geometry, bluntness and wall condition, and on the frequency and the orientations of the 
disturbances. However, when these conditions are known the computations of the boundary layers and the 
stability characteristics can be performed without much difficulty. 

The first process is called receptivity and it depends on several parameters including the flow features 
near the leading edge, the stability characteristics of the boundary layers and the types of external 
disturbances that impinge on the boundary layer. Since the phase speed of the acoustic disturbances 
synchronizes with the phase speed of the first modes near the lower branch of the neutral stability curve it 
is expected that acoustic disturbances will be very efficient in generating the instability waves in 
supersonic boundary layers. This observation led to numerous investigations on the interactions of 
acoustic waves with supersonic boundary layers. The investigations can be divided into two groups. One 
is the forced response of the parallel and non-parallel boundary layers to free stream acoustic 
disturbances 2 ' 3 ' 4 ' 5 . In these studies, it was found that acoustic waves excite disturbances inside the 
boundary layer that are much larger than those in the free stream. The computations also revealed that at 
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some critical Reynolds numbers the incoming acoustic waves are completely absorbed by the boundary 
layer. It was postulated that the first mode emerges downstream of this region with the initial amplitude 
given by the forced response calculations. The second group of investigations are concerned with the 
generation of the first and the second modes near the leading edge region due to the interaction with the 
free stream acoustic disturbances 6-7-8-9 . Fedorov and Khokhlov 6 investigated the excitation of the first and 
the second modes by the longitudinal acoustic waves using asymptotic theory. They showed that the 
growing boundary layer diffracts the incoming acoustic wave and generates eigen solutions inside the 
boundary layer. In a companion paper 7 they investigated the generation of the instability waves due to 
incoming acoustic waves at finite angles of incidence. It was identified that at small incidence angles the 
incident and the reflected waves are diffracted by the boundary layer resulting in the formation of the 
eigen solutions. At large incidence angles the scattered acoustic waves near the leading edge are 
diffracted by the boundary layer resulting in the generation of the eigen solutions. The receptivity of 
supersonic boundary layers is investigated numerically 8 by imposing a periodic Mach wave near the 
leading edge of the lower surface of a flat plate. It was shown that the first mode instability waves are 
generated in the upper surface of the flat plate and forced response was created in the lower side. 

At high Reynolds numbers the neutral point is located very close to the leading edge within one 
wavelength of the acoustic wave and the synchronization of the acoustic wave and the first mode occurs 
closer than this. Hence the flow features near the leading edge such as effects of bluntness; shocks and 
non-parallelism may become very important in the generation of the instability waves. These effects were 
not included in the aforementioned analysis or in the computations. In our previous studies 10- U- 12 , the 
generation and the evolution of three-dimensional disturbances induced by slow and fast acoustic 
disturbances and isolated roughness in a supersonic boundary layer over blunt flat plates and wedges at a 
free stream Mach number of 3.5 were numerically investigated by solving the full three-dimensional 
Navier-Stokes equations. It was found that instability waves are generated within one wavelength of the 
acoustic wave from the leading edge. The wavenumber first decreases from the acoustic value and then 
slowly increases to the instability wave value, which is smaller than the acoustic wavenumber. The 
computed receptivity coefficients based on the pressure fluctuations at the wall near the neutral points are 
about 1.20 and 0.07 for the slow and the fast acoustic waves respectively. 

Most of the receptivity analysis and computations are performed for flows over flat plates. Since the mean 
flow is uniform in the spanwise direction one Fourier mode with constant spanwise wavelength is 
considered in the analysis. In flows over axi-symmetric geometries the azimuthal wavenumber has to be 
an integer due to the circular periodicity in the azimuthal direction. These azimuthal wavenumbers for the 
unstable waves in supersonic boundary layers are large in the order of 10 to 30. When these wavenumbers 
are converted to wavelengths by dividing the circumferential length by the wavenumber, the wavelengths 
become very small near the nose and increase gradually in the axial direction. It was found in an earlier 
investigation 10 of receptivity of supersonic boundary layers to acoustic disturbances over flat plates that 
the instability waves are generated near the leading edge. The question is then how do the acoustic 
disturbances in the free stream with long spanwise wavelengths generate these small wavelengths or the 
large wavenumber disturbances inside the boundary layer over axi-symmetric bodies. 

There were several relevant transition experiments performed at NASA Langley in the Mach 3.5 
Supersonic Low-Disturbance tunnel. Boundary-layer transition data on a flat plate and on a cone along 
with free stream noise levels and the power spectral distribution of the free stream noise are presented in 
Ref. 13. As a first step, the computations are performed for the same conditions as in the experiment. The 
objectives are to understand the receptivity process near the leading edge of an axi-symmetric cone and to 
estimate the receptivity coefficient of the instability waves generated near the leading edge. A schematic 
diagram of the computational set up is depicted in Fig. 1 . 
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Governing Equations. 

The equations solved are the three-dimensional unsteady compressible Navier-Stokes equations in 
conservation form in cylindrical coordinates 
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Here (x, r, 9) are the cylindrical coordinates, ( u , v ,w) are the velocity components, p is the density, and p 
is the pressure. E is the total energy given by 

u 2 + v 2 + w 2 

E = e + , 

2 

e = c v T , p = pRT. (3) 

Here e is the internal energy and T is the temperature. The shear stresses are given by 
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The viscosity (p) is computed using Sutherland’s law and the coefficient of conductivity ( k ) is given 
in terms of the Prandtl number Pr. The variables p, p, T and velocity are non-dimensionalized by their 
corresponding reference variables p„, p„ T x and ^RT, respectively. The reference value for length is 
computed by Jvx 0 / U x , where x 0 is a reference location. For the computation, the equations are 
transformed from the physical coordinate system (x, r, 6) to the computational curvilinear coordinate 
system ,rj , £) in a conservative manner. The corresponding Cartesian coordinates are (x, y, z) with z 

and y-axis oriented along 6 = 0 and 90 degree lines respectively, Fig. (1). 


Solution Algorithm 

The governing equations are solved using a 5th order accurate weighted essentially non-oscillatory 
(WENO) scheme for space discretization and using a third order, total variation diminishing (TVD) 
Runge-Kutta scheme for time integration. These methods are suitable in flows with discontinuities or high 
gradient regions. The governing equations are solved discretely in a uniform structured computational 
domain where flow properties are known point wise at the grid nodes. In a given direction, the spatial 
derivatives are approximated to a higher order at the nodes, using the neighboring nodal values in that 
direction. The resulting equations are then integrated in time to get the point values as a function of time. 
Since the spatial derivatives are independent of the coordinate directions, multi dimensions can be easily 
added to the method. It is well known that approximating a discontinuous function by a higher order (two 
or more) polynomial generally introduces oscillatory behavior near the discontinuity, and this oscillation 
increases with the order of the approximation. The essentially non-oscillatory (ENO) and the 
improvement of this, WENO methods, were developed to maintain the higher order approximations in the 
smooth regions and to eliminate or suppress the oscillatory behavior near the discontinuities. These are 
achieved by systematically adopting or selecting the stencils based on the smoothness of the function that 
is being approximated. Reference 14 explains the WENO and the TVD methods and the formulas. 
Reference 15 gives the application of the ENO method to the N-S equations. Reference 16 describes in 
detail the solution method implemented in this computation. 

At the outflow boundary, characteristic boundary conditions are used. At the wall, viscous conditions 
are used for the velocities and a constant temperature condition is employed for the temperature. The 
density at the wall is computed from the continuity equation. In the azimuthal direction, symmetric and 
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periodic conditions are used at the boundaries. In the mean flow computations, the free-stream values are 
prescribed at the upper boundary that lies outside the bow shock. In the unsteady computations, the 
acoustic perturbations are superimposed On the uniform mean flow at the upper boundary. The procedure 
is to first compute the steady mean flow by performing unsteady computations using a variable time step 
until the maximum residual reaches a small value ~10' n . A CFL number of 0.15 is used near the nose 
region and 0.50 is used downstream. The next step is to introduce unsteady disturbances at the upper 
boundary of the computational domain and to perform time accurate computations to investigate the 
interaction and evolution of these disturbances downstream. There exists a singularity along the symmetry 
axis r = 0. This singularity is removed by taking the limit of the governing equations at r = 0. 

The symmetric acoustic field that impinges on the outer boundary is taken to be in the following 
form. 


r f ~ ia,„.x±iB nr z+i£ I 

Pac = R eal{p a(; e ~ } 

(7) 

+Real{ p gce to-**v~z-i*«y-t<* } . 

Here a ac , /3 ac , e ac are the acoustic wavenumbers, and a> is the frequency of the acoustic disturbance. The 
wavenumber in the y-direction e ac determines the incident angle of the acoustic waves and in this paper 
computations are performed for zero incident angle, e ac = 0.0. 

Results 

Computations are performed for supersonic flows over a 5-degree half-angle cone with blunt leading 
edges. Table 1 gives the flow parameters and Fig. 1 shows the schematic diagram of the computational 
set up. The nose region of the wedge is modeled as a circle of the form 

(x-r 0 ) 2 + y 2 = r 0 2 . (8) 

Here r 0 is the radius of the leading edge bluntness. The circular nose is smoothly merged with the cone by 
a tenth order polynomial. Unsteady simulations are performed for nose radii r 0 = 0.001 and 0.01 inches at 
a unit Reynolds number of 1.0*10 5 6 /inch. Stability computations also are performed for nose radius r 0 = 
0.001, 0.01 and 0.05 inches to investigate the effects of bluntness on the instability of supersonic 
boundary layers over axi-symmetric bodies. These parameters give the Reynolds numbers based on the 
nose radius equal to 10 3 , 10 4 and 5*10 4 . 


Table 1 Flow parameters for the wind tunnel model 


Free stream Mach number: M x = 3.5 

Free stream Reynolds number: Re ot =1.0*10 6 /in. 

Free stream density: p x = 2.249* 10" 2 lbm/ft 3 

Free stream pressure: p x = 1 87.74 lbf/ft 2 

Free stream velocity: U x = 2145.89 ft/s 

Free stream temperature: T X =\56A2 R 

Free stream kinematic viscosity: =1.7882*10" 4 ft 2 /s 
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Wall temperature: = adiabatic condition 
Prandtl number: Pr= 0.72 
Ratio of specific heats: y= 1.4 


Length scale 


V V- 


v = 5.892* 10 ~ 5 ft. (x 0 = 0.5 in.) 


The boundary layer thickness at x 0 = 0.5 in.: 8 0 = .0055 in. 
Non-dimensional frequency F=l*10" 5 is equivalent to 41.0 kHz 

The non-dimensional frequency F is defined as r = — — , 


where /is the frequency in Hertz. 


The grid is generated using analytical formulae. The grid stretches in the i] direction close to the wall 
and is uniform outside the boundary layer. In the direction, the grid is symmetric about the leading edge 
and very fine near the nose and is uniform in the flat region. The grid is uniform in the azimuthal 
direction. The outer boundary outside of the shock follows a circle near the nose region with its vertex 
located a short distance upstream of the nose and follows a parabola downstream of the nose to capture 
the boundary layer accurately. The computational domain extends from x = -0.015 to 10.0 inches in the 
axial direction. Calculations were performed using a grid size (2001*251*91) and (3001*351*136) 
depending on the wavelength of the acoustic disturbances. Due to the very fine grid requirement near the 
nose, the allowable time step is very small and the computations become very expensive to simulate the 
unsteady computations in the entire domain at once. To overcome this, calculations are performed in two 
steps. First, the computations are done near the nose region with a very small time step. Second, the flow 
properties in the middle of this domain are fed as inflow conditions for the second larger domain and the 
computations are carried out with a larger time step. 


Linear instability 

The local linear parallel stability analysis is performed for the axi-symmetric similarity boundary layers 
over the 5 -degree half angle cone at a free stream Mach number of 3.5. The conditions downstream of the 
shock on the surface of the cone are: 


M s =3.355 

T s = 1 .06 T,„ 

p s = 1.15 

p s = 1.16 p„ 

Re s /in =1.075*10 6 


Here, M s , T s , p s and p s are the inviscid surface Mach number, temperature, density, and the pressure 
respectively. Two forms of disturbances are used to perform the stability computations. 
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<Kr,x,z,t) = $(r)e ias+lpz - l " n 
(p(r,x,0,t) = ~4>(r)e ia5+,me - ,mt 

In the first form, (3 represents the spanwise wavenumber and in the second form m, which is an integer, 
represents the azimuthal wavenumber. The wavenumbers (3 and m are related by the relation 

m 
R 

where R is the radius of the surface. Hence if one assumes a constant /?, the azimuthal wavenumber m 
will vary along the axial direction and if m is assumed to be a constant the spanwise wavenumber (3 will 
vary along the axis. In axi-symmetric boundary layers the mathematically correct form for the disturbance 
is the second equation. However, if one neglects the curvature effects in the local stability analysis both 
forms will yield the same stability results. As the equation shows, near the nose region where the radius 
of the cone approaches zero the wavenumber /3 becomes very large or the spanwise wavelength becomes 
very small. These observations are very important when one investigates the receptivity of instability 
waves due to external disturbances in axi-symmetric boundary layers. How do the external acoustic 
disturbances with large spanwise wavelengths generate these very small wavelength instability waves 
near the leading edge? The possibilities are that the spanwise acoustic waves are modulated near the sharp 
nose region and the instability waves are generated in this region or the acoustic waves interact with small 
roughness near the nose region and the instability waves are generated. In this paper we consider the first 
possibility and investigate the interaction of acoustic waves with a cone in supersonic boundary layers. 


The linear stability results and the N-Factor computations for the similarity boundary layer over the 5 -deg 
cone are shown in figures 2 and 3. Figure 2 (a) shows the neutral stability diagram for different wave 

angles 0, 45, 60 and 70 degrees in the ^F,-^Re t j plane. Figures 2(b) and (c) show the N-Factor diagrams 

for the most amplified frequencies. Figure 2(b) shows the results for constant non-dimensional spanwise 
wave numbers [> and figure 2(c) shows the results for constant azimuthal wave numbers m. The most 
amplified frequency is about F = 2.0-2.5*10" 5 and the spanwise wave number of the most amplified wave 
is about /3 = 0.05. The azimuthal wave number of the most amplified wave is about m = 20. In the 
experiment the transition occurred at a Reynolds number of 3000 and this yields an N-Factor of 10 for a 
constant /3 and a factor of 8.5 for a constant m. The critical Reynolds number is about 700 and the critical 
Reynolds number for the most amplified disturbances is about 1000. To understand the length scales 
involved it is better to provide the wavenumbers in dimensional quantities. Figures 3(a) and (b) show the 
variation of the axial and spanwise wavelengths in inches and the growth rates (-cq) for the 

most amplified wave with the frequency F = 2.0* 10' 5 . Figure 3(a) shows the results for a constant 
spanwise wave number of (5 = 0.05 and figure 3(b) shows the results for constant azimuthal wave 
numbers of m = 10 and 20. The non-dimensional spanwise wavenumber (5 = 0.05 corresponds to a 
dimensional wavelength of 0.089 inches. Figure 3a shows that for a constant /3 = 0.05 the wavenumber m 
varies from 3 at a Reynolds number of 700 to about 50 at a Reynolds number of 3000. The axial 
wavelength X x increases from 0.168 inches to 0.210 inches and the waveangle increases from 60 to 65 
degrees. Figure 3(b) shows that for a constant azimuthal wavenumber of m = 10 the spanwise wavelength 
increases from 0.032 inches to 0.46 inches and the axial wavelength first decreases from 0.20 inches to 
0.183 inches and then increases to 0.212 inches. The waveangle decreases from 80 to 20 degrees. For a 
constant spanwise wave number of m = 20 the spanwise wavelength increases from 0.016 inches to 0.23 
inches and the axial wavelength first decreases from 0.23 inches to 0.192 inches and then increases to 
0.207 inches. The waveangle decreases from 85 to 40 degrees. 
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Acoustic waves 

The governing equations and the solutions for plane acoustic waves are given in Ref. 10. The acoustic 
pressure p is in the form 


P = 


~ iia^x+e^y + B 7-cot) 

v e ac at 

r ac 


The dispersion relation among the acoustic wavenumbers a ac , ft ac , £ ac and the frequency o> is given by 

(a ac U 0 -aj) 2 = (aj + ftj + £ ac 2 )al 
For zero incident angle e ac = 0, the wavenumber a ac can be expressed as 


OOCOSd_ 


a = 


(U 0 cosd z ±a 0 ) 

Here 9 = tan -1 ] — —I is the sweep angle and the plus sign corresponds to the fast moving wave and the 

K c) 

n 

minus sign corresponds to the slow moving wave. The corresponding phase speeds are C = U 0 


cos 


The wavenumber for the fast moving wave is a ac < 


co 


U 0 + a 0 


and for the slow moving wave 


a > 


00 


U 0 a 0 


and the sweep angle is limited by 0_ < cos | — |. At M = 3.5 the sweep angle of the slow 


moving wave is limited to 73.39 degrees. 


Table 2 Values of « ... and wavelength for different spanwise wavenumbers 


ft ac 

dac 

Kc, (in.) 

^ac,x (in*) 

X x (in.)m=10 

Neutral 

X x (in.)m=20 

Neutral 

0.00 

0.0198 


0.2244 



0.01 

0.0207 

0.4440 

0.2145 



0.025 

0.0274 

0.1776 

0.1847 



0.05 

0.0309 

0.0888 

0.1436 



0.10 

0.0455 

0.0444 

0.0976 



0.15 

0.0603 

0.0296 

0.0736 







0.1951 

0.2133 


Table 2 gives the non-dimensional wavenumbers a ac , ft,,, and the dimensional wavelengths T ac>x , >.. 1C Z for 
the slow acoustic waves. The axial wavelengths for the neutral disturbances for the azimuthal 
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wavenumbers m = 10 and 20 are also given in the table. It is seen that the shortest possible spanwise 
acoustic wavelength is about 0.0296 inches and the corresponding axial wavelength is about 0.0736 
inches. The wavelengths of the neutral disturbances are about 0.20 inches, which is about three times 
longer than the acoustic wavelength. Hence a conversion of both the spanwise wavelength and the axial 
wavelength has to occur to excite the instability waves inside the boundary layer. 


Mean flow and transition for the Cone 

Figure 4 shows the mean flow density contours computed using the WENO code. The figures 4(a-d) 
show the results for the cone at different nose radii r 0 = 0.001, 0.01 and 0.05 inches. The computations are 
performed at a unit Reynolds number of 1.0*10 6 /inch which yield Reynolds numbers based on the nose 
radius of 10 3 , 10 4 and 5*10 4 . Figure 4(a) shows the density contours in larger domain while figures 4(b) 
and (c) show the flow field near the nose region. The leading edge shocks are located approximately at 
0.0002, 0.002 in. and 0.02 inches upstream of the leading edge. The density profiles at different axial 
locations are plotted in Figs. 5(a)-(d) for the different bluntness cases r 0 =0.001, 0.01 inches in the 
physical and similarity coordinates. The compressible Blasius similarity profile is also included for 
comparison. The difference between the density profiles with the bluntness and the similarity profiles 
increases with the bluntness. For the bluntness case of r 0 =0.001 inches the boundary layer profiles 
approach the similarity profiles at about x = 2.0 inches while for the larger radius of 0.01 inches it 
approaches the similarity profiles at about 6.0 inches. 

Figure 6 shows the growth rate and the N-Factors for the most amplified disturbances computed using 
the profiles obtained from the numerical solution for different bluntness r 0 = 0.001, 0.01, 0.05 inches. For 
comparison, the stability results for the Blasius similarity profiles, which model a sharp leading edge, are 
also shown in the figure. The frequency and the azimuthal wave number for the most amplified wave are 
about F = 1.75*10' 5 and m = 25 for the similarity profiles. For the smaller bluntness r 0 = 0.001 and 0.01 
the N-Factor curves and the growth rate curves remain closer to the similarity curve. For the larger 
bluntness case r 0 = 0.05 the growth rates become smaller and the N-Factor curves move further to the 
right. The most amplified frequencies and the azimuthal wave numbers are ( 1 .75* 10" 5 , 25), ( 1 .75* 10' 5 , 25) 
(0.60*10‘ 5 , 25) for r 0 = 0.001, 0.01 and 0.05 inches respectively. This shows that the frequencies of the 
most amplified disturbances become smaller with increasing bluntness. The growth rate curves remain 
similar to the Blasius profile for all three bluntness cases. Boundary-layer transition data on a flat plate 
and on a cone, and free stream noise levels and the power spectral distribution of the free stream noise are 
presented in Ref. 9. The data shows the transition Reynolds number for a sharp cone in quiet conditions is 
about 9*10 6 . If this is used in the correlation to obtain the N-Factor at transition, then Fig. 6 gives an N- 
factor of 9.0 for a constant azimuthal wave number. The transition Reynolds numbers obtained using the 
N-factor of 9.0 for different bluntness cases are summarized in Table 3 and plotted in Fig. 7. The 
transition Reynolds number for blunt flat plates and blunt wedges were presented in Ref. 12 and are also 
included in the Fig. 7 and in the Table 3 for comparison. The ratio between the transition Reynolds 
number with bluntness and the transition Reynolds number for the similarity profile, (Re T ) b /(Re x ) similarity is 
about 1.05, 1.13, 2.10 respectively for Re,. = 10 3 , 10 4 , 5*10 4 . For the flat plate these factors are about 1.05, 

I. 30, 2.48, 3.49 respectively for Re b = 10 2 , 10 3 , 5*10 3 , 10 4 and for the wedge these factors are 1.42, 7.01, 

II. 96 respectively for Re b = 10 3 , 10 4 , 4*10 4 . This implies that the stabilizing effects of bluntness are much 
stronger in flows over flat plates and wedges than in flows over cones. Previous experiments 17- 18 and the 
stability calculations 19 showed that the transition Reynolds number for a blunt cone at a Mach number of 
8 with nose Reynolds numbers of 30,000 increased by a factor of 1.7~2.0 compared to a sharp cone. This 
agrees with the calculations shown in the Table 3. 
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Table 3 Transition Reynolds number for the blunt flat plates, wedges and cones 


Nose 
radius 
r 0 in. 

Reynolds 
number 
based on 
nose 
radius 

Transition location X 
(in.) 

Transition Reynolds 
number *10 6 

Ratio 

(Re x ) b / (Re T ) Sim j| arity 

Flat 

plate 

Wedge 

Cone 

Flat 

Plate 

Wedge 

Cone 

Flat 

Plate 

Wedge 

Cone 

0.0 

0 

11.00 

10.40 

9.37 

11.00 

10.40 

9.37 

1.00 

1.00 

1.00 

0.001 

10 3 

14.12 

14.38 

9.80 

14.12 

14.38 

9.80 

1.30 

1.42 

1.05 

0.01 

10 4 

38.37 

73.00 

10.60 

38.37 

73.00 

10.60 

3.49 

7.02 

1.13 

0.01 (4 
mil) 

4*10 4 


28.50 



114.00 



10.96 


0.05 

5*10 4 



19.60 



19.60 



2.10 


Generation of instability waves from a harmonic source 


The computer code was modified from the Cartesian coordinates to cylindrical coordinates. To validate 
and to investigate the generation and the evolution of instability waves in axi-symmetric boundary layers 
a simulation with a harmonic point source located on the surface is performed. A blowing/suction 
distribution in the form 


v = Cnax sin 


~ X start ) 

( X end ~ X start) / 


cos/770 sin(atf) 


is prescribed for the radial velocity component across a narrow slot on the surface. Here x smm x end are the 
starting and ending locations of the slot, m is the azimuthal wavenumber, at is the frequency and v max is 
the maximum velocity non-dimensionalized by the free stream velocity. In the present computations x start 
= 0.70 in., x end = 0.80 in., m = 10, F=2.00*10' s and v max = 0.00001. The boundary layer thickness at the 
axial location x = 0.75 inches is about 0.007 inches. Figure 8 shows the evolution of the unsteady density 
fluctuations generated by the harmonic point source at a fixed time. Figure 9a shows the instantaneous 
pressure fluctuations along the wall and Fig. 9b shows the amplitude of the pressure fluctuations along the 
wall in a log scale. Figure 9b also includes the results from the parabolized stability equations (PSE) 
computations obtained for the same mean boundary layer profiles. The figures clearly show the initial 
generation and the eventual exponential growth of the instability waves inside the boundary layer. The 
growth of the disturbances agrees very well with the PSE results about ten boundary layer thicknesses 
downstream of the source. Following the PSE results up to the location of the source, the initial amplitude 
of the instability waves near the source can be estimated. From these values the receptivity coefficients 
defined by the initial amplitude of the pressure fluctuations at the wall near the harmonic source non- 
dimensionalised by the maximum source velocity can be evaluated. 

sP wall ) n 


c 


recpt,p waU 


( 15 ) 


The computed receptivity coefficient is 


C , =0.50 

recpt.p waU 
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Figure 10 shows the eigenfunction distributions for the density, p, and the azimuthal velocity, w, 
obtained from the simulation and the PSE computation at the axial location x = 2.0 inches. The agreement 
between the simulation and the PSE is excellent. 


Interaction of a three-dimensional acoustic wave with the axi-symmetric boundary layer. 

After the mean flow is computed, three-dimensional plane acoustic disturbances are introduced at the 
outer boundaries. The simulations are performed for a non-dimensional frequency of F=2.00*10' 5 and for 
spanwise wavenumbers j5 ac = 0.05 and 0.15. The amplitude of the acoustic wave is p ac l p x = 1.0 * 10 -6 . 
The computations are performed for zero incident angles. Figures 11 (a) and (b) show the results for the 
evolution of the unsteady density fluctuations obtained from the simulation for the slow wave /3 ac = 0.15 
at a fixed time in the (x, z) and (x, y) planes respectively. It is to be noted that the azimuthal angle 0 = 0 
degrees lies in the (x, z.) plane and 0 = 90 degrees lies in the (x, y) plane. The constant phase lines in the 
(x, y) or 6 = 90 degree plane is perpendicular to the cone axis. Figure 12 depicts the unsteady density 
fluctuations in the cross plane (z, y) at an axial location x = 0.5 inches. It is interesting to see the 
scattering of the acoustic field by the cone and the different ways the flow field is perturbed along the 0 = 
0 and 90 degrees lines. The constant phase lines are parallel to the surface in the 0 = 0 direction and the 
acoustic field basically glances the cone without much distortion. In contrast, the constant phase lines are 
diffracted by the cone in the 0 = 90 direction and the acoustic field is very much distorted. These 
observations will reflect the way the boundary layers are perturbed by the three-dimensional plane 
acoustic waves. It is to be noted that the acoustic disturbances in the outer part of the domain in fig. 11 
are decaying due to lack of resolution in this region due to the expansion of the domain. 

Figures 13 (a) and (b) show the plan view of the density fluctuations on the surface of the cone in (x, rO) 
planes. Figure 13 (a) depicts the flow field between the azimuthal angles (-90, 90) degrees so 0 = 0 line is 
in the middle of the region and figure 13(b) depicts the field between (0, 180) degrees so 0 = 90 line is in 
the middle. As it is seen, the disturbances evolve differently along the 0 = 0 and 90 degrees lines. Figures 
14(a)-(d) show the pressure fluctuations on the surface along different azimuthal directions 0 = 0, 42, 70 
and 90 degrees. Along the 0 = 0 line, the wavelength remains the same as the acoustic field up to x = 1.5 
inches and it gradually starts to increase downstream. However, along the 0 = 90 degree line the 
wavelength increases starting from the leading edge and the disturbance amplitude starts to increase 
beyond x =1.5 inches. It is also seen that the wavelength modulation and the amplitude growth spreads 
in the azimuthal direction as the disturbances evolve downstream. 


Figure 15 shows the pressure fluctuations on the surface along the azimuthal direction at a fixed time at 
different axial stations x = 0.10, 0.25, 0.50. 1.24 and 2.00 inches. The figure shows that near the nose, x = 
0.1 inches, the variation in the azimuthal direction is similar to a simple cosine wave and near x = 0.25 
inches small wavelength disturbances start to appear in the middle region 0= 90. Further downstream, the 
wavelengths or the number of waves m increases and these small wavelength disturbances spread towards 
0 = 0 and 180 degrees. To investigate how the disturbances are evolving, a Fourier analysis is performed 
in the azimuthal direction. The disturbances are decomposed into the following form. 


p(x,r,6,t ) 


2) Pn,( X ^) e ' me [<?■'"* +C.C 


m=-M 
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Figure 16 shows the amplitude of the Fourier component \p m (x,r)\ as a function of the azimuthal 
wavenumber m for the pressure distribution on the surface of the cone at different axial locations x = 
0.042, 0.27, 0.52, 1.00, 1.49, 2.72 inches. As was discussed earlier, the spectrum is narrow and is 
concentrated near m = 0 for small x values and it widens as the axial distance increases. The amplitudes 
of the disturbances with azimuthal wavenumbers in the range of 8 to 12 start to increase beyond x = 1.0 
inches. Figure 17 shows the variation of the amplitude of the spectrum |p m (x,r)| along the axial 
direction x for different m = 0, 10 and 20. The figure shows that the amplitude of the disturbances with m 
= 0 decreases and the amplitude for m = 10 increases sharply up to x = 0.5 inches and plateau up to x = 
1.0 inches and increases exponentially downstream. However, the disturbance with high azimuthal 
wavenumber m =20 increases first and then decreases with increasing axial distance. The linear stability 
analysis predicted that the azimuthal wavenumbers for the unstable waves range between 5 to 30. It is not 
clear why the disturbances with large m values are not amplified but the disturbances with m values in the 
range of 6 to 14 are amplified. Figure 17 also includes the results from the parabolized stability equations 
(PSE) computations obtained for the same mean boundary layer profiles. The figure clearly shows the 
eventual exponential growth of the instability waves inside the boundary layer. However, the initial 
amplitude of the instability waves is very small compared to the free stream acoustic disturbances. The 
initial amplitude from the PSE computations at x = 0.62 inches is about 5.0* 10' 9 compared to the free 
stream acoustic disturbance amplitude of 1.0* 10' 6 . Hence, the receptivity coefficient is very small value 
5.0* 10" 3 . 

Similar to the previous case, figure 18 show the plan view of the density fluctuations on the surface of the 
cone in (x, r&) planes for the case /3 ac = 0.05. Figure 18 (a) depicts the flow field between the azimuthal 
angles (-90, 90) degrees and figure 18(b) depicts the field between (0, 180) degrees. The observations are 
similar to the previous case with the disturbances evolve differently along the 6 = 0 and 90 degrees lines. 
Figures 19(a) and (b) show the pressure fluctuations on the surface along the different azimuthal 
directions 6 = 0 and 90 degrees. Surprisingly the disturbances did not grow along the 6 = 90 degree as 
was observed in the previous case with large /3 ac = 0.15. These questions need further investigation. 


Discussion and Conclusions 

The boundary layer receptivity process due to the interaction of three-dimensional slow and fast 
acoustic disturbances with a blunt 5 -degree cone and the evolution of disturbances inside the boundary 
layer are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 
10 6 /inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations 
using the 5 th -order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization 
and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. 

The results show that bluntness has a less stabilizing effect on the stability of the axi-symmetric 
boundary layers compared to two-dimensional boundary layers. The transition Reynolds number for the 
cone increases slowly up to a nose Reynolds number of 10000 and increases sharply at higher nose 
Reynolds numbers. This suggests that up to a nose Reynolds number of 10000, the cone could be 
assumed as a sharp cone for the stability analysis. The transition Reynolds number for a flat plate with 
leading edge thickness of 0.01 inches is about 3.5 times larger than that for the Blasius boundary layer, it 
is about 7.2 times for the wedge and it is about 1.13 times for the cone. Even at a very high nose Reynolds 
number of 50,000 this ratio is 2.1 for the cone. 

The unsteady simulations revealed that the interaction of a three-dimensional plane acoustic wave with a 
cone perturbed the flow field inside the boundary layer differently depending on the azimuthal direction. 
In regions where the constant acoustic phase lines are parallel to the surface the disturbances have the 
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same wavelength as the acoustic waves. In the regions where the constant phase lines are perpendicular to 
the surface the boundary layer is highly perturbed. The axial wavelengths of the disturbances increase and 
disturbances with high azimuthal wavenumbers appear in this region. Further downstream, these 
disturbances spread in the azimuthal direction. The amplitude of the excited disturbances is very small 
compared to the free stream acoustic disturbances. The receptivity coefficient is about 5.0* 10' 3 . The 
simulation with high spanwise acoustic wavenumber showed that these high azimuthal wavenumber 
disturbances generated near the nose region become the instability waves and grow exponentially 
downstream. However, the simulation performed with the small spanwise acoustic wavenumber did not 
show any instability waves even though the modulation of the disturbances near the nose region appear 
similar. Another point of concern is that according to the linear stability analysis the disturbances with 
azimuthal wavenumbers in the range of 5 to 30 are unstable at the Reynolds number investigated in this 
paper. However, the simulation showed that the disturbances with high azimuthal wavenumbers of 20 did 
not amplify. These question will be investigated further. 
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p ac = a exp( i (a ac x±(5 ac z+e ac y-0jt))+c.c 



Figure 1. Schematic diagram of the computational model. 


(a) 




(c) 



Figure 2. Neutral stability (a) and N-Factor diagrams (b), (c) for the boundary layer over a 5-deg 

cone. M„ = 3.5. 
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(a) Constant |3=0.05 , F= 2.0*10' 5 (b) Constant m=10, 20, F = 2.0*10' 5 




Figure 3.Variation of the axial and spanwise wave lengths, the growth rates and the azimuthal wave 
numbers for the most amplified wave F = 2.0*10' s . (a) constant p=0.05, (b) constant m =10, 20. 





Figure 4. Contours of the density for flow over a 5-degree cone with a blunted leading edge at 

M=3.5 
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Figure 5. Mean density profiles at different X locations. 
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Figure 6. N-Factors and the growth rates for the most amplified disturbances for different 
bluntness r 0 = 0.001, 0.01, 0.05 in. at M = 3.5. 



Figure 7. Transition Reynolds number for the flat plate and the wedge for different bluntness b 

0.0001, 0.001, 0.005, 0.01 in. at M = 3.5. 
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Figure 8. Contours of the unsteady density fluctuations due to a harmonic point source at x = 0.75 

inches. E=2.0*10' 4 , m= 10. 




Figure 9. Amplitude of the pressure fluctuation at the wall and comparison with the PSE. 

E=2.0*10' 4 , m = 10, r 0 = 0.001 inch. 



Figure 10. Eigenfunctions for the density and azimuthal velocity fluctuations and comparison with 
the PSE calculations at X= 2.0 inches. E=2.0*10' 4 , m = 10, r 0 = 0.001 inch. 
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(a) 9= Odeg., in the (x, z) plane 


(b) 9= 90 deg., in Ihe (x, y) plane 
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Figure 11. Contours of the unsteady density fluctuations due to the interaction of a 3D slow 
acoustic wave with a 5-deg cone with a blunted leading edge. F = 2.00*10' s , /? = 0.15, incident 
angle 0.0. 
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Figure 12. Contours of the unsteady density fluctuations in the y_z plane at x = 0.5 in. due to the 
interaction of a 3D slow acoustic wave with a 5-deg cone with a blunted leading edge. F = 2.00*10' s , 
yS = 0.15. 
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Figure 13. Contours of the unsteady density fluctuations due to the interaction of a 3D slow 
acoustic wave with a 5-deg cone with a blunted leading edge in the (x, r0) plane. F = 2.00*10' s , /? 


0.15. 
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(a) (b) 






Figure 14. Pressure fluctuations on the wall along different azimuthal angles. r 0 = 0.01 in. F = 
2.00*10 5 , yS = 0.15. 



Figure 15. Pressure fluctuations on the wall at a fixed time along the azimuthal direction at 
different axial locations. r 0 = 0.01 in., F = 2.00*10' s , /3 = 0.15. 


22 of 24 

American Institute of Aeronautics and Astronautics 


37 th A1AA Fluid Dynamics Conference and Exhibit, June 25-28, 2007, Miami, Florida 



Figure 16. Spectral distribution at different axial locations. F = 2.00*10'% yS = 0.15. 



Figure 17. Variation of the amplitude of the spectrum along the axial direction for different m. F 

2.00*10' s , /5 = 0.15. 
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Figure 18. Contours of the unsteady density fluctuations due to the interaction of a 3D slow 
acoustic wave with a 5-deg cone with a blunted leading edge in the (x, r0) plane. F = 2.00*10' s , /? 
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Figure 19. Pressure fluctuations on the wall along different azimuthal angles r g = 0.01 in. F = 
2.00*10' 5 , yS = 0.05. 
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